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Abstract. We report results of molecular-dynamics simulations of a model polymer melt consisting of 
short non-entangled chains in the supercooled state above the critical temperature Tc of mode-coupling 
theory (MCT). To analyse the dynamics of the system we computed the incoherent, the collective chain 
and the collective melt intermediate scattering functions as well as their space Fourier transforms, the van 
^ ■ Hove correlation functions. In this first part of the paper we focus on the dynamics in the /3-relaxation 

' regime. The final structural relaxation, the a-relaxation, will be studied in the following second part. The 

results can be summarized as follows: Without using any fit procedure we find evidence for the space-time 
factorization theorem of MCT in real and reciprocal space, and also for polymer specific quantities, the 
Rouse modes. The critical amplitudes in real space are determined directly from the simulation data of 
the van Hove correlation functions. They allow to identify the typical length scales of the /3-dynamics, 
and illustrate that it is a localized process. In a quantitative analysis the wave-vector dependences of 
the /3-coefficients, i.e., of the non-ergodicity parameter, the critical amplitude, and the next-to-leading 
, order correction coefficients, are studied for all correlators. The /3-coefficients show indications of polymer 

■ specific effects on the length scale of the chain's radius of gyration. The agreement between simulation 

and the leading-order MCT description is found to be good in the central /3-regime. Next-to-leading order 
corrections extend the validity of the MCT approximations to a greater time window and become more 
important at large wave-vectors. 

PACS. 64.70.Pf Glass transitions - 61.25.Hq Macromolecular and polymer solutions; polymer melts; 
swelling - 61.20.Ja Computer simulation of liquid structure 

O 1 Introduction to and below Tc additional relaxation mechanisms grad- 
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ually become more important than the cage effect. The 
A distinguishing property of glass forming liquids is the approximate inclusion of these processes in the theory (ex- 
strong increase of the structural relaxation time, which tended MCT) avoids the absolute freezing at Tc The 
precedes the glass transition, i.e., the solidification of the glass former remains liquid down to Tg, where it falls out 
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' liquid on the experimental time scale at the glass tran- of equilibrium. According to MCT, Tc is a characteristic 

' sition temperature Tg ||l|,^,||. Understanding the micro- (thermodynamic) temperature of glass forming hquids. 

scopic origin of this remarkable slowing down of the dy- mi r ^ ^ ^ i-^-^ 

^. ^ , „ , 1 1, • 1, 1 he existence of an additional characteristic tempera- 

namics represents one ot the most cnallengmg problems ,. ,ji ^ ^ x^ii 

, , , . . , ^. ° ture IS suggested by recent attempts to connect the struc- 

m condensed matter physics. A solution to this problem , , , °. -^^ ^ n ^ ^ ■ ^ j i -^.i 

, , , 1 • /n T^rn\ iHl n n Ri tural glass transition to mean-neld spm-glass models with 

was proposed by mode-coupling theory (MCI I P.q,q,u1. j- j i. / rFl F 



The theory predicts the existence of a critical temperature 



discontinuous order parameter (see for a review). 

„ „ , „ ; 1- 1 .1 1 • ^^. 1- ^ This approach is motivated by two observations: First, the 

ic, well above ig, at which the dynamics qualitatively , , . , . r j i ■ j -u j u 

, ' , i , i m 1 ■ • 1 • long-time dynamics ot these models is described by the 

changes. Above, but close to ic the dynamics is domi- • i t i i t j_- r ^ ^- i i 

^ ? , ^1111- r 1 1 -1 idealized mode-coupling equations oi schematic models, 

nated by the mutual blocking oi a particle and its neigh- r i -u-^ j • ^ -i- ^ ^ 

, rr o, \ mi - iT 1- , /r \ i heretorc, tucy cxhi bit a dvuamic trausitiou at a tcmpera- 

bours ( cage effect I. ihis blocking creates (free) energy , i ■ i! • -j .-c j -^i ^ o j ^ j 

, . ^ , . , /, , • n -7 111. ture which IS identified with Jc- Second, they also undergo 

barriers which would become mfanite and lead to a cease ^ ^- , -i- ^ • i i ^ . . 

r , , ,1 .r.i rr;i; -la static transitiou to a spm-glass phase at a temperature 

of any structural relaxation it the cage effect determined m i i m t ^i ■ ^ i m m m i 

1 . 1 /• 1 1- 1 Tvir^mN TT 1 J^s below ic- In the interval Jg < J < Jc the dynamics 
the dynamics alone (idealized MCi I. However, very close •,, - ji ^-ui i^ rii 
' I J jg determined by an exponentially large number ot local 
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a long time before escaping to another one by activated 
processes. Below Tg the number of minima becomes finite 
and the system is trapped in one of them. 

At present, it is not clear whether this analogy be- 
tween spin and structural glasses complements or chal- 
lenges mode-coupling theory. A recent application to ex- 
periments jlTf shows that the theory overestimates Tg. It is 
placed in the liquid phase close to Tc- On the other hand, 
this application as well as other experimental or simula- 
tion studies provide evidence for the existence 
of a critical temperature Tc, as predicted by MCT. Such 
evidence has motivated extensions of MCT to treat ori- 
entational degrees of freedom ir2i|T^ , p^ , vibrations |p^ , 
[l^ or non-equilibrium systems |17f|. Furthermore, many 
computer simulations |l8|l have been undertaken to un- 
derstand better the dynamics above and below Tc (dy- 
namic heterogeneity ||l|,|^,|l[|2|,|2|,||, connection be- 
tween the potential energy landscape and structural relax- 
ation [||,|6[|2^,||,|9|, physical aging [|o[|3l|), but also to 
further test the theory on other systems than simple liq- 
uids. Some examples are diatomic molecules p2|, p3|p^ ,|35[, 
ortho-terphenyl [^|3^, water |^,^,40,4l| or polymers 

With this and the subsequent paper |45[] (hereafter 
called part I and part II, respectively) we want to con- 
tribute to this research by extending our previous analy- 
sis of a simple model for a supercooled polymer melt. This 
analysis considered the relaxation behavior of the incoher- 
ent intermediate scattering function [ ^Of , the dynamics of 
the model under isobaric and isochoric conditions 
the interplay between the cage and polymer-specific ef- 
and correlated motion of the monomers close 
51|. The present work attempts to complement 
these studies by discussing results for the intermediate 
coherent scattering and van Hove correlation functions of 
both the polymer and the melt. Whereas part II deals 
with the final structural relaxation of the a-regime, the 
present part I describes the dynamics of the model in the 
/3-relaxation regime defined by MCT close to Tc. It is or- 
ganized as follows: Section || briefly introduces the model. 
Section || compiles the quantities studied and the theoret- 
ical background for the analysis. The following section ^ 
summarizes the results and the final section |^ contains our 
conclusions. 



2 Model 



= 1.1225(7 being the minimmn of [/l,t- Throughout this 
paper all quantities are measured in Lennard- Jones units: 
temperature and distances were measured in units of e/fce 

and (T, respectively, and time in units of (mcr^e)^^^, with 
mass m set to unity. 

Along the polymer backbone an additional FENE (fi- 
nitely extensible nonlinear elastic) potential ||5^ was used 
to introduce bonds between adjacent monomers. 



k 



Ro 



(2) 



with Rq = 1.5 and k = 30. The superposition of the LJ 
and FENE potentials leads to an effective bond potential 
with a sharp minimum at rbond = 0.9606, making bond 
crossings impossible. 

These model parameters yield two incompatible length 
scales: rbond = 0.9609 and r,„iii — 1.1225. "Incompatible" 
means that it is not possible to arrange the beads in a 
regular fee (or bcc) lattice structure if one requires that 
the beads take precisely these distances which yield the 
pairwise minima of the potentials. Another feature of the 
model is that the chains do not become stiffer with de- 
creasing temperature. This is reflected by the nearly con- 
stant end-to-end distance R^ — 12.3 ± 0.1 and radius of 
gyration i?g — 2.09 ±0.01 in the investigated temperature 
range . The flexibility along the backbone of the chain 
and the incompatibility of rbond and rmin prevent crys- 
tallization of the melt ||5^. The static structure exhibits 
the typical features of an amorphous material in the tem- 
perature range studied (see Fig. |^). This is an essential 
premise for investigating glassy dynamics. 

The simulations were done in two steps: First, the vol- 
ume of the simulation box was determined in a constant 
pressure simulation ai p — 1. Then, this volume was kept 
flxed and the simulations were continued in the canoni- 
cal ensemble using the Nose-Hoover thermostat. Periodic 
boundary conditions were applied in all three spatial di- 
rections of the cubic simulation box. A box typically con- 
tained between 110 and 120 chains with = 10 monomers 
each. We simulated at T = 0.46, 0.47, 0.48, 0.50, 0.52, 
0.55, 0.6, 0.65, and 0.7. While the lowest simulated tem- 
perature T = 0.46 is slightly above Tc ~ 0.45, the crit- 
ical temperature of the mode couphng theory 
the highest temperature T = 0.7 leaves the t emp erature 
range where the /3-dynamics can be observed [t46|] . 



In this section we briefly describe the model underlying 
our simulation. A more detailed description can be found 
in|7|.^ ^ 

We investigated a bead-spring model of linear polymer 
chains. All monomers interact via a truncated and shifted 
Lennard- Jones (LJ) potential given by 



C/Lj(r) := 



4e 



(7)" -(7) 



0, 



r < 2ri, 
r > 2rn 



(1) 



where the constant C = 0.00775 is chosen such that the 
potential vanishes continuously at r = 2rmin, rmin — 2^/^(t 



3 Simulated Quantities and Theoretical 
Background 

The static structure of a glass is almost indistinguish- 
able from that of a liquid. The two phases mainly dif- 
fer in the values of their relaxation times: The relaxation 
time in the liquid phase is small, but becomes unmea- 
surably large in the glassy phase. A means to investigate 
this crossover from fast to strongly protracted structural 
relaxation are dynamic correlation functions. In this pa- 
per we present our flndings for the coherent intermediate 
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scattering and the van Hove correlation functio ns i n the 
/3-relaxation regime. In the following subsection 3.1 these 
functions are defined and their calculation from simulation 
data is briefly described. 

In the vicinity of Tc mode-coupling theory (MCT) pre- 
dicts a two-step relaxation process for correlation func- 
tions which couple to density fluctuations So me 
MCT results are compiled in the second subsection \i.2\ 
as well as approximations for the single chain scattering 
function, which are motivated by polymer theory p4l. 



3.1 Analysed Quantities 

The coherent intermediate scattering function of the melt 
is defined by 



Ml 

S,{Q) 



(3) 



Here, Sq{0) = Sq is the collective static structure factor 
and Sq{t) is given by 



^«(^) ^= tf(E E • ['■^(*) - ^^(0)] } 



M M 



M\ 



(4) 



.1 j=i 



where ri{t) is the position of the ith monomer at time t 
in the melt, which contains M monomers in total, and (•) 
denotes the canonical ensemble average. 

Due to isotropy all discussed quantities depend on the 
modulus g = |q| of the wave vector q only. In the following 
we will use a superscript "s" for incoherent (self) and a 
superscript "p" for chain (polymer) quantities. Coherent 
quantities of the melt do not carry a superscript. 4>q{t) 
stands for any scattering function, where the superscript 
"x" indicates the different correlators. In the case of the 
chain scattering function, (j>^{t) is defined by equations (||) 
and (^ when using N, the number of monomers per poly- 
mer, instead of M and summing over all pairs of monomers 
belonging to the same chain. Setting i = j in equation (^ 
defines 4>q{t), the incoherent scattering function. 

Since the system is homogeneous and isotropic, one 
can simplify the formulas in a way that is more suitable 
for computation. For instance, Sq{t) may be calculated by 



^ M M 
Sq{t) = — ^^COs[qTi(t)]^cos[q 



M 



M 

sin[q • ri{t)] ^ sin[q • rj(0)] 



(5) 



i=l 



which allows computation with a number of trigonometric 
operations of order M, instead of order M^. For S^{t) and 
Sq{t) analogous expressions are obtained. 

The set of g-values for which the scattering functions 
were calculated is g e {1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 6.9, 7.15, 
8.0, 9.5, 11.0, 12.8, 14.0, 16.0, 19.0}. These values cover 
the maxima and minima of the static structure factor (cf. 



Figure Especially for small q there are only very few 
reciprocal vectors complying with the periodic boundary 
conditions. This leads to unsatisfactory statistics in the 
case of the scattering function of the melt (for S^{t) the 
statistics is M/N times better). In order to improve the 
statistics the scattering functions were averaged over an 
interval of width A around q. So, whenever quoting a func- 
tion value <pq{t), this actually means 



ci\qe[q-A/2,q+A/2]} Vcj 

\{ci\qe[q-A/2,q + A/2]]\ 



(6) 



where the sum runs over all allowed reciprocal vectors. A 
was chosen such that it is as small as possible, but provides 
sufficient statistics, and it was checked that the averaging 
did not introduce significant deviations. 

Insight in the dynamics in real space can be obtained 
by analysing the van Hove correlation functions, which are 
the Fourier transforms of the respective scattering func- 
tions. For the melt the van Hove correlator is defined by 
(see e.g. ii) 



M M 



(7) 



i=i j=i 



G(r, t) is proportional to the probability density for find- 
ing a particle in the volume element dr. It is useful to 
separate G{r,t) into self and distinct parts, G(r,t) = 
Gs(r,0 + Gd(r,t) with 



G'^(^'^) = ]^(E'^[^ + ^''(0)-^^W])' (8) 
1=1 

M M 

- -(5]5]<5[r + r,(0) -r,(t)]) . (9) 



i=i j=i 
Hence, we have at t = 

Gs(r,0)=<5(r), Gd(r,0) = ^^.g(r) , (10) 
where we introduced the pair correlation function g{r) 

A I M 



9{r) 



V 



M{M - 1) 



^^(5[r + r,-r,;] 



(11) 



i=i j=i 



The denominator M(M — 1) is the total number of mo- 
nomer pairs in the melt. For polymers all quantities are 
defined in the same way, with M replaced by N, and indi- 
cated by the superscript "p" . Again, we introduce a short 
hand notation, Gx(r, i), to denote self and distinct van 
Hove functions. In the case of isotropic systems G(r, t) 
depends on r — |r| only. After discretization the van Hove 
functions can be computed from histograms of pair dis- 
tances in the melt, taking periodic boundary conditions 
into account. 
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3.2 Theoretical Background 



This section is split into two parts. The first part 3.2.1 



summarizes those MCT predictions which are relevant for 
the subsequent analysis of the scattering and van H ove 
correlation functions, whereas the second part 3.2.2 dis- 



cusses approximations, suggested by the theory of poly- 
mer dynamics , for the coherent scattering function of 
a chain. 



3.2.1 Mode-Coupling Theory 



t f), and so to times on the scale to-, where |0^(t) — 

/ri « 1- 

The time-evolution of (j)q{t) in the /3- regime can be 
written as [p6l 



ux-i-a 
fa 



9(t) 



2b 



(16) 



"'q'-O 



-2a 



The subsequent discussion mainly follows references pq , 
|57[ , where one can also find quantitative results for a sys- 
tem of ideal hard spheres. 

The central prediction of MCT is the existence of a 
critical temperature Tc in the temperature regime of the 
supercooled liquid. In the vicinity of asymptotic results 
can be derived for small values of the separation parameter 



Tr-T 



T, 



{12) 



The constant C is system dependent and typically of order 
1 |5^. The following equations are predictions in leading 
order of \a\. 

In the idealized MCT the separation parameter de- 
termines the temperature dependence of the two relevant 
time scales of the theory, i.e., of the /^-relaxation time t^, 



tc 



in 



|CT|l/2a ' 



and of the universal a-timescale f, 



T>Tr. 



(13) 



(14) 



In equations (|T^) and (|lj) represents a matching-time 
to the microscopic transient, and the exponents a (0 < a < 
1/2) and 7 are related to each other by 7 = l/(2a)-|-l/(26). 
Here, h denotes the von Schweidler exponent which is in 
turn connected to the critical exponent a by the exponent 
parameter A via 



r{i-af _ r{i + bf 



r{l-2a) r{i + 2b) 



The exponent parameter, and thus also a, b, and 7, de- 
pend on the glass former under consideration, but not on 
temperature or on the correlator 0^(t) studied. 

Close to Tc, (j>q{t) is predicted to decay in two steps: 
A first step leads to a plateau value, the so-called non- 
ergodicity parameter /™, and a second step leads off of 
it. The second step represents the initial part of the final 
structural a-relaxation, whose temperature dependence is 
determined by f . In this scenario the /^-relaxation regime 
corresponds to the time-window between to and f {to <C 



where t :~ t/t„. The first two terms of equation (|T^ con- 
stitute the so-called "factorization theorem" . The name 
derives from the property that the function g(i) carries 
the whole temperature and time dependences, while the 
dependence on space and on the correlator enters only via 
fg'^ and the critical amphtude h^. Therefore, g{t) is called 
(3 master function. Its shape is solely determined by the 
exponent parameter A. 

The factorization theorem is the MCT result in lead- 

1/2 

ing order of |cr| (it is of order \a\ ' ). The second and third 
lines of equation (16) represent corrections (of order \a\) to 
it. The third term of the sum ( |l6| ) is the next-to- leading or- 
der long-time correction, which becomes important in the 
late /3-regime, whereas the last term gives the correspond- 
ing correction for the early /3-relaxation. These corrections 
violate the factorization property of the leading order re- 
sult due the g-dependence of K^. Expressions for Kg, 
and n have been derived in ||5^,^. For the analysis of the 
simulation data, it is only important that these quantities 
are given in terms of the static structure at Tc and can 
thus be considered as independent of temperature (in the 
same way as f^^ and h^). 

To perform this analysis we rewrite equation ( p^ ) by 
combining those theoretical quantities which cannot be 
determined separately in a fit procedure. Since the micro- 
scopic time to is unknown, only the product 



,x fit 



(17) 



can be fitted. Furthermore, we introduce long- and short- 
time correction coefficients, and B^, defined by 



Bq - 



(18) 
(19) 



(15) This yields the following expression for equation (g^ 



fit 



-9ii) 



fit 

^A^gP" 

f2a q 



A^B^t-^'' 



(20) 



In this equation only t„ varies with temperature according 
to equation (13). All other quantities are (asymptotically) 
independent of temperature. Futhermore, note that A^ 
and Bg satisfy the following linear relationship 



(21) 
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Both correction coefficients, and should thus ex- 
hibit the same g-dependence. Equation ( |2l| ) can either be 
used as a test of this prediction or for calculating the 
short time coefficient from the long-time coefficients 
A^, which were easier to obtain from the simulation data 



(see section 4.2.5). 

From a previous analysis of the incoherent scattering 
function the exponents a, b and 7, the /3-relaxation 
times ta, the parameter B, and the f3 master function 



g(t) are known. In section 4.2.5 we want to exploit this 
knowledge as much as possible to extend the quantitative 
description of the /3-dynamics of (t)q{t) to coherent scat- 
tering. 



3.2.2 Coherent Chain Scattering Function and Rouse Model 

In a dense polymer melt excluded volume and hydrody- 
namic interactions are screened [|4|. If such long-range 
interactions are absent, the Rouse model is commonly as- 
sumed to provide a viable approximation for the dynamics 
of non-entangled chains [ ^8[|59|| . Although the entangle- 
ment length iVo has not been determined for the present 
model yet, extensive simulation studies of a closely related 
model suggest Nc w 32 [E^,^. Since our simulations are 
done with = 10 (< A^cX^the Rouse model should apply. 
Therefore, we want to compare the coherent chain scat- 
tering function with the predictions of the Rouse theory 
in the following. 

A basic assumption of the Rouse model is that every 
monomer experiences a local random force which is Gaus- 
sian distributed. Since the displacement ri{t) — rj(0) is a 
linear function of these random forces, it is also a Gaus- 
sian random variable 1 54 . This implies that the canonical 
average of the exponential in equation can be written 
as 



(exp{iq. [r,(i)-r,(0)]}) 

= exp|-^([r,(i)-r,(0)]')| (22) 



so that 



^ N N . 2 ^ 

1=1 7=1 ^ ' 



(23) 



Note that the time-dependence of this "Gaussian appro- 
ximation" for S'P(i) is completely given in terms of the 
mean-square displacements ^i{t) — rj(0)]^) and that the 
g-dependence is quite different from that of the factoriza- 
tion theorem (|l|). Both S\{t) and ([r.,(t) - rj(0)]2) can 
be determined independently in the simulation. By com- 
paring the results we are able to test to what extent the 
displacements can be considered as Gaussian distributed 
random variables. 

A further step consists in expressing (|rt(t) — rj(0)]^) 
by the Rouse modes Xp(t), defined by Q, 



1 ^ 

Xp(t) := — ^r,(t)cos 



p7r(i - 1/2) 



TV 



(24) 



where p (= 0, 1, . . . , A^ — 1) denotes the mode index. For 
our systems the Rouse modes were analysed in This 
analysis showed that the Rouse-mode correlation func- 
tions, (Xp(t) • Xg(0)), are (to a very good approximation) 
orthogonal at all times, i.e., (Xp(t) • Xq(0)) oc (5pg. When 
using this property the displacements may be expressed 
as 



33(0+4^(^^(0)) 



[r,(t)-r,(0)] 

337r(i- 1/2) 



X I cos 



N 



pnU - 1/2) 



A^ 



(25) 



X cos 



>7r(i - 1/2)' 




p^(j-l/2)" 


N 


cos 





where 33 (i) denotes the mean-square displacement of the 
chains' centre of mass, and 



(Xp(0 . Xp(0)) 



.N-1) 



(26) 



(Xp(0) • Xp(0)) 
are the normalized time-correlators of the Rouse modes 



Again, both sides of equation (|25[) can be obtained inde- 
pendently from the simulation. We call the Gaussian ap- 
proximation ( p3| ) with displacements calculated via equa- 
tion (^) "Rouse approximation" . 



4 Simulation Results 
4.1 Statics 

For the interpretation of the intermediate scattering func- 
tions it is important to know the static structure factors 
Sq and S^. Sv, was computed in [^6j and is compared with 
S'P in Figure Both Sq and exhibit an amorphous halo 
at (7max ~ 7 corresponding to the nearest neighbour shell. 
As q grows, Sq oscillates around 1, approaching 1 for the 
highest accessible q. These features are characteristic of 
an amorphous structure and present at all temperatures. 
With decreasing temperature the maxima and minima be- 
come sharper, and the position of the first maxima shifts 
to slightly larger q, as the melt becomes denser. 

At small q the structure factor of a chain, SJj, can be 
fairly well described by the Debye function p4| 



•^DcbycC?) 



2N 



q^R 



-9 -K, 



+ q^Rl - 1 



TV 



1 + 



(27) 
(28) 



Equation (^8|) represents an approximation to the De- 
bye function (p7|), which is accurate to about 15% [p4| . 
The Debye function assumes a Gaussian distribution of 
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Fig. 1. Static structure factors of the melt Sq and of the 
chain 5, . Sq is shown at three temperatures covering the range 
of investigated temperatures Tc ~ 0.45 < 0.46 < T < 0.7. 
Sq depends only weakly on temperature and exhibits typical 
features of a liquid structure at all T. Sq remains unchanged in 
this temperature region and is shown together with the Debye 
approximation (p7|). Reciprocal lengths corresponding to the 
end-to-end distance Rc — 3.5 and the radius of gyration _Rg ~ 
1.45 are also indicated. 



the monomer-monomer distances, an assumption which 
is certainly not justified at higher q probing smaller dis- 
tances where the exact form of the potentials matters. 
Therefore, deviations between simulation and theory are 
expected for q > 2n/Rg. At g « 7.5 the simulated has 
a maximum caused by the nearest neighbour shell. Com- 
pared to Sq it is shifted to larger q because the distance 
between bonded monomers in a polymer is shorter than 
between non-bonded nearest neighbours. varies only 
very slightly in the studied temperature range. Thus, the 
data at T = 0.46 are a representative example. 

The corresponding quantities in real space, the pair 
correlation functions of the melt and the chain, g{r) and 
g^{r), are shown in Figure |^. The main features remain un- 
changed in the whole temperature range, but at lower tem- 
perature maxima and minima become more pronounced. 
The overall appearance resembles closely the pair corre- 
lation function of hard spheres and simple Lennard-Jones 
hquids (see e.g. |]55|). 

However, there are also characteristic differences. In 
a simple liquid there is only a single nearest neighbour 
peak, whereas the polymer model exhibits two peaks: one 
at Tbond = 0.9609 and another one at rmin = 1-1225 
corresponding to the two length scales of the model 

The peak at Tmin reflects the sharp minimum of the 
bond potential. The minimum of the Lennard- Jones po- 
tential is less steep. Hence, one expects a wider peak. 
Only monomers in the same chain which are not direct 
neighbours can contribute to the peak at rmin for g^{r). 
Therefore, it is much lower than that of g{r) where sur- 
rounding polymers can also contribute. The subsequent 
next-nearest shells cause the periodic sequence of maxima 
and minima in g{r), whose amplitude gradually decreases 



300 



^200 

3. 




Fig. 2. Pair-correlation functions of the chain (top) and of 
the melt (bottom) for the lowest temperature T = 0.46 and 
the highest temperature T = 0.7. The inset in the upper 
panel enlarges the region of r > 1.1. The vertical lines are at 
^bond = 0.9609, the preferred bond length, at rmin = 1.1225, 
the preferred distance of non-bonded monomers, and at r = 
iig = 1.45, the radius of gyration. The circles in the top panel 
represent the approximation ( |29[ ) for 5Dcbyc('')- 



as r grows. There is also a weak, broad second neighbour- 
shell maximum in g^{r) around r ~ 1.8. It is situated at 
smaller r than the corresponding maximum of g{r) at r = 
1.95, presumably because next neighbours are closer in the 
chain than in the melt, roughly by r^am — ''bond — 0.1616. 
The inset in the plot for g^ shows a magnification of the 
comparison with the Fourier-transform of the approxima- 
tion (^8|) for the Debye function {y is the volume of the 
simulation box) 



3Dcbyc('') 



27V 



V 1 



(27r)3/2(7V- 1) Rl r 



exp 



V2r 



(29) 



Equation ( p9|) describes gP(r') quite well for r > 3, where 
details of the monomer-monomer interactions do not mat- 
ter. 



4.2 Dynamics 
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Fig. 3. Coherent scattering function 05 (t) at T = 0.65 
wave vector moduli q = 1, 2, 3, 4, 5 
Symbols show the simulation results 
the Gaussian approximation (Eq. (Es 



for 

6.9, 9.5, 14, and 19. 
Dashed lines represent 
), solid lines the Rouse 



approximation, i.e., Eq. (M) together with Eqs. (Ea) and (Eq) 
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Fig. 4. Same plot as in Figure ^, but at T = 0.48 and for q = 1, 
2, 3, 4, 5, 6.9, 9.5, and 14. The simulation results are repre- 
sented by the symbols, whereas the dashed and solid lines show 
the Gaussian (Eq. (0)) and Rouse approximations (Eq. ( p3| ) 
with Eqs. (Ea) and p6[)), respectively. 



4.2.1 Intra-chain Coherent Scattering Function and Rouse 
Analysis 

In this section we compa re the Gaussian and Rouse ap- 
proximations (see Section 3.2.2| ) with the simulation data 
for the intra-chain coherent scattering function at different 
q. For T = 0.65 Figure ^ shows that both approximations 
agree with one another and describe the correlators quite 
well. Deviations occur for q 4 and 5 around t — 0.5 and 
for all especially for q > 6.9, in the a-regime. These 
observations suggest the following conclusions: The agree- 
ment between the Gaussian and the Rouse approximations 
shows that equation (|25| ) represents an accurate descrip- 



tion of ([ri(i) — rj(0)]^). This means that the Rouse modes 
are, to a very good approximation, orthogonal at all times 
for our model at high temperatures. Thus, the difference 
between simulation and theory must be attributed to non- 
Gaussian distributed displacements ri{t) — rj{Q). The in- 



fluence of this non-Gaussian character depends on q and 
t, but in general makes the simulated scattering functions 
decay more slowly than the theory. The same conclusion 
was also drawn from atomistic simulations of polyethylene 
[ p2| and polybutadiene 63, 6^. 

At T = 0.48, both the Gaussian and Rouse approxi- 
mations poorly describe the a-relaxation for q > 3, and 
fail completely for large q already at the beginning of 
the /3-regime (see Figure Q). Although a two-step relax- 
ation behaviour is reproduced, the collective effects in the 
polymer are not properly taken into account. Since the 
Gaussian approximation fails, it can be concluded that the 
forces acting on the monomers do not give rise to Gaus- 
sian distributed displacements as temperature approaches 
Tc. On the other hand, the figures show that both the 
Gaussian and Rouse approximations are almost indistin- 
guishable at all q. Thus, the Rouse-mode formula ( |2^ ) 
for the displacements remains accurate even in the super- 
cooled regime close to Tc (this was also checked by compar- 
ing the displacements from simulation data directly with 
equation (|2F 



4.2.2 Dynamics of the melt in the /3-relaxation regime in 
g-space 

We want to test the validity of the factorization theorem 
(cf. equation (|l6|)) 



by computing the function |^ 



fa 



9{t/ta) 



(30) 



g{t/t„) - g{t'/t„) ^ ^^^^ 



<p-{t') g{t"lt„)~g{t'/t,) 



(31) 



If the factorization property holds, R^{t) depends only on 
the fixed times t' and t" and on temperature (via t^), but 
not on the correlator "x" or on q. The times t' and t" 
can be chosen arbitrarily in the /^-regime. For numerical 
stability it is, however, advisable to take values at the 
beginning and the end of the plateau region in order to 
obtain a large denominator. An advantage of equation (|3l| ) 
is that it represents a simple test which works directly with 
the simulation data. No intricate fit procedure is involved. 
This approach has therefore been pursued in reciprocal 
and in real space in several other studies |6^,^,^ (see 
also [ |68||69|] f or comparable experimental tests). 

Figure ^ shows R^{t) computed from the scattering 
functions at T = 0.46 with t" = 0.610 and t' = 86.43. This 
temperature is already so close to Tc (~ 0.45) that devia- 
tions from the ideal MGT prediction in the a- regime, i.e., 
from the power law behaviour (|lj) of the a-timescale, oc- 
cur |6H4|,||J|L Such deviations are also found in other 
simulations [ p2| , |3^ ,p^ or in experiments, see e.g. and 
are usually attributed to additional relaxation channels 
which are not taken into account by idealized MGT, but 
become dominant as temperature decreases towards and 



8 



M. Aichele, J. Baschnagel: /3-Dynamics of Glassy Polymer Melts 



below Tc Nonetheless, Figure 5| illustrates that there 
is an intermediate time- window of about 2.5 decades where 
the data for incoherent and coherent scattering collapse 
onto a g-independent master curve, while they splay out 
at both short and late times. We tested that the master 
curve is the same for all correlators "x". Very similar re- 
sults are obtained for higher temperatures, with times ap- 
propriately chosen in the plateau-region for each temper- 
ature, as long as a two-step relaxation is observed (i.e., for 
T < 0.52 iQ). The observation that the factorization the- 
orem is satisfied for temperatures where ergodicity restor- 
ing processes already violate the validity of equation ( p^ 
suggests that it remains valid also for T <Tc (as expected 
theoretically 

Furthermore, Figure ^ qualitatively confirms predic- 
tions for the higher order corrections to the factorization 
theorem: According to equation (|2l|), the short- and long- 
time correction coefficients have the same g-dependence. 
These corrections become important before and after the 
central /3-regime. This means that the top curve before the 
collapse onto the master curve is also the top curve after 
the collapse, the second from top before t" is the second 
from top after t' , and so on. This behaviour is illustrated 
with dashed and dot-dashed lines in Figure || for some 
q. Except for the coherent scattering functions, which ex- 
hibits oscillations at small q in the early /3-regime (see 
part II of this paper for a more detailed discussion of this 
point ) the prediction is seen to be fulfilled. One finds 
that the correction coefficients for the incoherent scatter- 
ing function, and B^, grow monotonously with g, in 
agreement with calculations for an ideal hard sphere sys- 
tem |5^. Similar findings were reported in ]6q| for the 
incoherent scattering function of a binary Lennard-Jones 
fluid. 



4.2.3 Dynamics of the Rouse modes in the /3-relaxation 
regime 
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Fig. 5. Test of the factorization theorem for all scattering 
functions at T = 0.46 by plotting Rq{t), defined by equa- 
tion (^l|), vs. t. From top to bottom the three panels show 
the results for incoherent, coherent-chain and coherent-melt 
scattering. In the /3-regime all curves collapse onto a sin- 
gle master curve. By definition, Rq{t" = 0.610) = 1 and 
i?J(t' = 86.43) = 0. 



Motivated by Figure H it is tempting to apply equation ( pT| ) 
also to polymer specific quantities, such as the Rouse mo- 
des. Therefore, we define i?p°"^°(i) in analogy to Rg{t) 
with 0^(t) replaced by <Pppit) (see Eq. (H)). Figure | 
shows R^^^'^{t) (for all Rouse modes p) and compares 
it to R^it) at T = 0.48. Note the different scale for the 
y-axis, causing the spread of Rq{t) to appear larger than 
in Figure |[ Again, t" and t' are chosen at the beginning 
and the end of the plateau. 

As with the scattering functions, all Rouse correlators 
collapse onto a master curve in the /3-plateau region. Al- 
though the splaying out at short and long times is much 
weaker than for the incoherent scattering function, the fig- 
ure suggests that the factorization property also holds for 
the Rouse modes. However, the master curve of R^°^^'^{t) 
lies slightly above that of Rg{t). Hence, one can speculate 
that an extension of MCT to polymer melts might yield 
a f3 master function for the Rouse modes, which is close, 
but perhaps not identical to g{t). 



On the other hand, this conjecture has to be considered 
with some reservations, since the Rouse modes exhibit a 
two-step relaxation with a rather high plateau value (the 
non-ergodicity parameter is presumably larger than 0.9 for 
all p) 14 9|] . This proviso is suggested by a recent compara- 
tive analysis of depolarized-light scattering, dielectric-loss 
and incoherent neutron scattering results for propylene 
carbonate [|6j . The study shows that the data from these 
different techniques can be consistently described by a 
schematic MCT model which goes beyond asymptotic pre- 
dictions and incorporates corrections resulting from hop- 
ping processes and from the crossover to vibrations and 
to the a-process. By comparing this description with an 
analysis using only the asymptotic laws it was found that 
the leading-order results can be strongly masked for quan- 
tities with large non-ergodicity parameters. 

Finally, Figure ^ allows another observation. It reveals 
a (further) difference between i?^°""°(i) and i?^(i). The 
corrections to the factorization property for the Rouse 
modes behave oppositely to those of R^ (t) : The order from 
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Fig. 6. Rl{t) and i?^°"'="=m for T = 0.48. is defined 

analogously to equation ( |3l| ) by replacing <p^{t) by the Rouse 
modes $pp{t). The times t' and t" axe chosen as t" — 0.988 and 
t' = 21.97. The order from top to bottom of the K,{t) curves 
outside the collapse region is conserved (see Figure H) , whereas 
the order of the curves is reversed (exemplified for 

p — 1 drawn as a thick dashed line). 

top to bottom of the curves before t" is reversed 

after t' for all R^°'^^'^{t) (only shown for p = 1), whereas 
it is maintained for the scattering functions (see also Fig- 
ure H) . 
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4.2.4 Dynamics of the melt in the /^-relaxation regime in 
real space 

The factorization property can as well be investigated in 
real space for the van Hove functions, for which one ob- 
tains by Fourier transformation of equation (pO|) 



Gxir,t) = Fx(r) + [t^/t^Y Hx{r)g{i) 



We now define 



i?x(r, t) 



Gx(r',<)-Gx(r',t') Hx{r') 



(32) 
(33) 



which was proposed in |£6|^ . We can choose any time t' 
in the /3-regime and an arbitrary value for r' . A possibility 
to fix r' is to take the position of the nearest neighbour 
maximum at t = for the distinct correlators and that of 
the maximum of 47rr^Gs(r, t) in the centre of the /3-regime 
for the self cor relator: r'f = 0.9575, = 1.1025, and 
Tg = 0.13253. Since the correlators decay faster at the peak 
positions than at other r, the denominator of equation ( p3[ ) 
is largest at r^. This is again favourable for numerical 
stability. Furthermore, Rx{r'x,t) = 1 by definition. If the 
factorization property is satisfied, the results for Rx{r,t) 
should fall onto a master curve at any time t in the window 
of the /3-relaxation. The master curve is specific for each 
correlator, but independent of T, as the critical amplitude 
Hx{f) does not depend on temperature. 



Fig. 7. Test of the factorization theorem for the van Hove cor- 
relators at T = 0.48. Rx{r,t) is depicted for 9 times from the 
interval t" = 0.988 < t < t' = 21.97 (same values for t" and 
t' as in Figure |^. For the distinct correlators the respective 
pair correlation functions (rescaled) are shown as dotted lines 
for comparison, r* is the zero of Rs{r,t). For 0.46 <T< 0.52 
all Rx{r,t) fall on a temperature independent master curve for 
each correlator. This master curve is the same for all temper- 
atures. In order not to overload the figure only the results for 
T — 0.48 are shown. The dashed lines correspond to t closest to 
t' where one can see numerical instabilities (note that Rx{r,t) 
is undetermined for t = t'). 



We calculated i?x {r, t) from the van Hove functions at 
various T. As an example. Figure shows the results for 
T = 0.48, which is the centre of the temperature range 
wher e a q uantitative /3-analysis was performed (see sec- 
4.2.5| and |46 ). As expected from the factorization 



tion 

theorem there is a characterisitic, temperature indepen- 
dent master curve for each correlator. The distinct cor- 
relators Hd{r) and H^ir) roughly follow the respective 
pair correlation functions if r > 1 (= monomer diame- 
ter). Comparable result s are also obtained in reciprocal 
space (see section 4.2.5| ). The spatial dependence of and 
the differences between different correlators can however 
be seen more clearly in real than in g-space (compare Fig- 
ures and H). The behaviour of Rd{r,t) for our system 
is qualitatively similar to that of a binary Lennard- Jones 
mixture (67[0 and of hard spheres . It is important 
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to note that the master curves of Figure |^ are only ob- 
tained for times t" <t < t' chosen from the intermediate 
time- window of the plateau (see Figures [s] and ^ . If one 
leaves this time-window towards microscopic times or to- 
wards late times of the a-process, the superposition of the 
data is no longer possible. The same observation was also 
made for the binary Lennard- Jones mixture | |6^ . 

Due to the factorization property the spatial variation 
of Hx{r) / Hx{r') provides information about the length 
scales which are involved in the /3- relaxation. Figure il- 
lustrates that all Hx{r) quickly vanish on the scale of a 
few interparticle distances. The most long-ranged decay is 
found for the critical amplitude of the melt's distinct cor- 
relator Hd{r) whose oscillations are damped out if r > 4. 
Thus, the /^-relaxation involves monomer rearrangements 
up to about the forth nearest neighbour shell. The dom- 
inant contribution comes from cooperative displacements 
of a monomer and its first and second neighbours. This 
local character of the spatial variation of Hx{r)/ Hx{r') is 
an illustration of the "cage effect" 

Let us consider the spatial variation of Hx{r)/Hx{r') 
in some more detail. The critical amplitude of the self 
correlator Hs{r) has a root at r* = 0.2323. This value is 
(almost) identical to \/6''sc — 0.2327, i.e., to the root of 
the approximation, 



H,{r') ~ (1 - r'2/6r2J exp ( - r'^/irl) 



(34) 



where r^c = 0.095 ± 0.005 |46|. Equation (g) is obtained 
by Fourier transformation of a Gaussian approximation for 
the critical amphtude /i^ (see Eq. (31b) of |Q). This ap- 
proximation provides a reasonable description of the sim- 
ulation data for ft,^ if q < (/max [0 • Although it cannot be 
quantitatively precise when r < 27r/(7max ~ 0.9, it qualita- 
tively reproduces the features of Hs{r) for r > r* (root at 
•\/6rsc, minimum around VTOrgc — 0.3) and proposes that 
the relevant length scale for the monomer motion in the /?- 
regime is Vgc- Numerically, the magnitude of rgc (= 0.095) 
is comparable to the Lindcmann melting criterion, stat- 
ing that a (crystalline) solid melts if a particle moves on 
average more than 0.1 of its own diameter around its equi- 
librium position. This suggests that the local motion of a 
monomer in the /3-regime resembles that of a particle in 
an amorphous solid which is about to melt. 

A similar conclusion was also drawn from a detailed 
analysis of the mean-first passage time in simulations of 
a binary LJ-mixture of small and large particles |2^. If 
temperature is close to Tc, the mean- first passage time 
of a large particle to reach a distance r from its origin 
grows most steeply for 0.21 < r < 0.25. Furthermore, 
the distribution of first passage times exhibits a power- 
law decay for this range of r (see Figs. 8 and 9 of pOt). 
These features are indicative of intermittency in particle 
motion for distances which are close to r* of the binary 
mixture (see Fig. 7b of |^) and of our model. Therefore, 
the analysis of Ref. pQ| supports the interpretation that 
the vanishing of iJs(r]at t* defines a length scale Vgc for 
transient particle localization in an amorphous structure, 
as also suggested in 0, for instance. 



Figure |^ shows that Hs{r) is positive if r < r*, but 
negative for r > r* . Due to equation ( |3^ ) the product 
Hs{r)g{t) determines Gs{r,t), i.e., the probability to find 
monomer displacements of size r in time t. Since the /3 
master function g(t) decreases monotonously p6| , we can 
write for the rate by which the van Hove functions decay 
with time 



dGx{r,t) 
dt 



= -{to/t^THx{r) 



dgii) 



dt 



oc ~Hx{r) . (35) 



Therefore, the probability for monomer displacements of 
size r in the /3-regime decreases most at the maximum po- 
sition of Hs{r) (r « 0.13) and increases most at its mini- 
mum position (r s:! 0.35), whereas it vanishes if r — )■ 1 (= 
monomer diameter). On the level of the individual mono- 
mer motion this should imply that a monomer is partially 
refiected back to its origin when attempting to penetrate 
into the zone initially excluded by its neighbours. This 
"caging" has been nicely demonstrated in a simulation of 
a (polydisperse) hard sphere system [ p2| by calculating 
the probability that a particle moves parallel to the di- 
rection of its preceeding displacement. It was found that 
this motion is on average oriented opposite to the initial 
direction if r < 0.8. 

The dynamics of the distinct correlators H^jr) and 
Hd{r) can also be interpreted by equation (35). Since 
H^{r) and Hd{r) are in phase with the respective pair- 
correlation functions if r > 1 (Figure equation ( |35| ) il- 
lustrates that the van Hove correlators G^{r, t) and Gd{r, t) 
change most at the positions where the probability to find 
another particle at t = was largest or lowest. The initial 
positions of high distinct-monomer density are strongly 
depleted, whereas the initial exclusion zones are popu- 
lated. So, other monomers also partially enter the region 
r < 1 that a monomer occupied at i = 0. In this re- 
gion the behavior of H^{r) and H^^r) resembles a mirror 
image of the shape of Hs{r) for r > r* . This is not un- 
reasonable. If the probability for a monomer displacement 
increases particularly at a distance r ^ 0.35 from its ini- 
tial position, the approach of other particles towards this 
position should be most pronounced around r « 0.7 due to 
the preferred intermonomer distances rbond — 0.9609 and 
''min — 1.1225, assuming a nearest neighbour replaces this 
particle. The comparison of H^{r) and H^ir) further sug- 
gests that bonded nearest neighbours penetrate the initial 
exclusion zone of a monomer more easily than non-bonded 
ones as H^{r) < H^^r) for r < 1. This should imply that 
adjacent monomers along the backbone of a chain tend 
to follow each other. On the other hand, the small values 
of both \Hl{r)\ and \Hd{r)\ for r < 0.7 indicate that the 
replacement of a monomer by its neighbours occurs only 
slowly with time in the /3-regime. The replacement process 
is slow because the monomer motion is in turn confined 
to displacements r < 1 by its neighbours. This is another 
evidence for transient "caging" . A microscopic analysis of 
this monomer motion is underway [^ll . 
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4.2.5 Quantitative /5-analysis 

In this section we turn to a detailed quantitative analysis 
of the /3-process. From previous work on the incoherent 
scattering function a lot of information is available 
which we want to exploit as much as possible to extend the 
description to coherent scattering. Our aim is to reduce 
any further fitting to a minimum. The starting point of 
the analysis is equation (pO|). In and in this work it 
was found that the asymptotic forrnulas for the description 
of the /3-relaxation are applicable for 0.47 < T < 0.52. 
Thus, this is the temperature range where the /3-analysis 
is carried out. 

The exponents a and b, and the constant B were de- 
termined in Ref. [EgI, 



0.352 ±0.010, 6 = 0.75 ±0.04, B 



0.476 ±0.060 . 

(36) 



Additionally, the P master function g{t) was numerically 
given, along with the /3-times for all temperatures. 
Thus, the complete time and temperature dependences 
in the /3-regime were known. 

Using the f3 master function the condition (?(tco) = 
allows us to determine the crossover time ico = tcota- ■ Now, 
we can read off the non-ergodicity parameters from the 
simulation data by 



(37) 



since higher order corrections are vanishingly small at tco 

Theoretically, /™ should not be temperature depen- 
dent, but we observed a slight dependence on T. So, we 
averaged over all temperatures to obtain the values shown 
in Figure |8[ f^'^ and /|"^ are very close to one another, be- 
ing different only at small q. For q > 11 all curves fall 
on top of each other. is in phase with Sq, as theoret- 
ically predicted for hard spheres p6| and also found in 
simulations of a binary Lennard- Jones liquid ||7|,0 and 
diatomic molecules |32] . Interpreted within the framework 
of MCT this means that modes at the next neighbour dis- 
tance (where Sq has its maximum ^max) decay very slowly 
and contribute strongly to the slowing down of the dy- 
namics of the melt. Interestingly, there is a weak shoulder 
at g « 4.5 corresponding approximately to the size of a 
chain, i.e., 27r/i?g « 4.35. A similar shoulder was also 
found in simulations of diatomic molecules , but is ab- 
sent for ideal hard spheres |5^. An explanation could be 
that the polymer coil may be viewed as a soft, diffuse par- 
ticle with radius i?g, which contributes on its length scale 
to the slowing down of structural relaxation. The same 
conclusion can be drawn from the a-relaxation times as 
discussed in part II This conjecture could be further 
tested by simulations with polymers of different length. 

In the vicinity of tco higher order corrections are neg- 
ligible and the critical amplitude (combined with a pref- 



actor to Eq. (|17|)) can be obtained by numerically 



evaluating 



' X fit 



dt[g{t)/t-] 




(38) 



Fig. 8. Non-ergodicity parameters f^"^, fg'^, and fg versus q. 
fg is in phase with the static structure factor, Sg, shown as a 
thin solid line. 



Again, the curves for /i^ and are very simi- 

lar, especially at large q as shown in Figure ^. The most 
prominent structure is observed for h^^. In the ideal hard 
sphere system there is a sharp minimum at Qmax, whereas 
hq exhibts only slight variations for larger q, which are 
roughly in phase with Sq . For diatomic molecules 
and water a distinct minimum was also found at g,nax- 
For q < gmax, hq first passes through a (more or less pro- 
nounced) maximum in these systems before it decreases 
(and eventually levels off) if q ^ 0. This is in contrast 
to our system where hq first shows a minimum and then 
grows for q — > 0. However, one has to be cautious with the 
results of coherent melt scattering at small q. On the one 
hand, the statistics is worse than at larger q, and on the 
other hand, we observed oscillations at the beginning of 
the /3-process in (j)q{t)m this g-range, which might be finite 
size effects (see Fig. and part II |^ for a more complete 
discussion). Despite this proviso it is interesting to note 
that hq exhibits a minimum at g = 4, which is close to 
the position where fq has a shoulder. So, one may specu- 
late that these features reflect the slowing down of a chain 
much in the same way as the shoulder and the minimum at 
<Zmax signal that of particles in a non-polymeric liquid. The 
minimum at gmax is visible here in form of a step at g « 7 
only. In contrast to the previously mentioned systems 

we find a clear maximum at q sa 8 > gmax- For 
q> 8, the curve oscillates with Sq (around /i^*^*), which is 
qualitatively in accord with results for hard spheres . 

Up to now, we have not considered higher order cor- 
rections to the factorization theorem, that is, the third 
and forth terms of equation ( pO| ) . They were found to be 
quite significant outside the centre of the /3-process [|6[ 
p7| . These higher order corrections cannot be read off di- 
rectly from the simulation data, but have to be deter- 
mined by fits. In order to extract them the values cal- 
culated from the factorization theorem (first two terms in 
Eq. (H^)) were subtracted from the simulation data. Next, 
we chose an interval around tco where fitting was per- 
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Fig. 9. Critical amplitudes fe^ and hf (Eq. (0)) 

togetiier with the static structure factor Sq (thin solid line) 
versus q. 
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formed. We took g-independent intervals beginning and 
ending approximately at the inflection points of the corre- 
lator before and after the plateau. The dependence on the 
chosen interval was not very strong, but noticeable. The 
short-time correction coefficients are in general hard 
to determine as the microscopic dynamics influences the 
results. The dependence of the early /3-relaxation on the 
microscopic dynamics of the simulation was demonstrated 
e.g. in p3. By comparing Newtonian and stochastic dy- 
namics it was found that the different short-time behavior 
does not affect the late-/3- and the a-process, but the early 
/3-process. The MCT-prediction for the approach towards 
the plateau was visible for stochastic microscopic dynam- 
ics only. Since the microscopic dynamics of our simulation 
is Newtonian, we cannot expect to obtain accurate results 
for from fits of the simulation data. 

Figure ^ shows that the long-time correction param- 
eters of the coherent melt scattering function Aq are (ap- 
proximately) in anti-phase relative to Sq. This is in qual- 
itative agreement with the findings for hard spheres |p6| 
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Fig. 11. Short-time correction coefficients Bq as a function of 
the long-time correction coefficients Aq (shown in Fig. |l^) for 
q = 1, 2, 3, 4, 5, 6, 6.9, 7.15, 8, 9.5, 11, 12.8, 14, 16, 19. MCT pre- 
dicts a linear dependence with slope — 4.4135 (Eq. (|2l|)), 
which is approximately satisfied at small q for the incoherent 
scattering function. For small q, extraction of Bq from simu- 
lation data is very difficult, and no smooth behaviour can be 
observed. At large q, all curves Bq{Aq) roughly agree with one 
another, but here a precise determination of Ag is hampered 
by the low value of the /3-plateau. 



and diatomic molecules 



. As in the case of the non- 
ergodicity parameters and the critical amplitudes, the val- 
ues of Aq and A^ are very similar, not showing much struc- 
ture. 

MCT predicts a linear relationship between A^ and 
B^ (cf. Eq. (^)). Figure [ill plots B^ versus A^ and com- 
pares the result with the theoretical prediction: a line with 
slope 1/B^ = 4.4135. Whilst the curve for B^(A^) exhibits 
linear behaviour with a slope close to 1/B^ as long as 
<Z < Qmax, the curve bends upwards at larger q. The coef- 
ficients of the coherent chain correlator behave similarly, 
except for very small q, where one observes deviations. 
Contrary to that, Bq{Aq) exhibits a very irregular struc- 
ture. There is no linear behaviour at all, presumably due 
to the problems mentioned above (oscillations at the be- 
ginning of the /3-regime for q < qmax, insufficient statistics) 
so that Bq cannot be fitted reliably. On the other hand, all 
curves almost coincide at large q, where the slope is sig- 
nificantly different from 1/B^. In this g-range, however, 
one has to take into account that the small plateau value 
(cf. Figure |^) renders the determination of the long-time 
correction coefficient Aq difhcult. This might lead to de- 
viations relative to the theoretical prediction. 

Now, we want to test whether the higher order correc- 
tions improve the description of the /3-process. Since our 
focus is on the coherent scattering functions, we show com- 
parisons of 4'qit) and the MCT approximations in leading 
and next-to- leading order {(j)^{t) is very similar to (j)q{t) 
which has already been discussed in Q|). In Figure |l^ 
(pq {t) is depicted for some q and for the highest and lowest 
temperature for which a quantitative /3-analysis was per- 
formed. As pointed out above, the short-time correction 
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Fig. 12. Coherent scattering function from simulation 

data (symbols) for q = 3.0, 6.9, 8.0, and 14.0 at T = 0.47 
(top) and T = 0.52 (bottom), in comparison with MCT fits 
in the /3-regime. The vertical lines are at t = fco (oc to-), the 
centre of the /3-regime. Lowest order approximations (factor- 
ization theorem) are drawn with dashed lines, solid lines are 
used for the next-to-leading order approximations with short 
time correction coefficients computed from long-time correc- 
tion coefficients (see Eq. (|2l|)). Note that different scales are 
used for the time axis in the upper and lower panels. While the 
factorization theorem gives a good approximation around tco, 
higher corrections must be considered at later times. As the fit 
results for fq shows a significant temperature dependence at 
higher q compared to small q, the temperature-averaged val- 
ues for fq used here lead to a visible off'set of the curves for 
q = 14.0. 



coefficients Bq could not be obtained reliably from the 
simulation data. Therefore, we calculated them from the 
more reliable long-time correction coefficients Aq, using 



-0.108 + 4.4135^„, where the slope is l/B^ and the 



offset is taken from the straight line in Figure [111. For both 
temperatures, the factorization theorem provides a good 
description if one is not too far away from tco- Especially 
in the late /3-regime next-to-leading order corrections im- 
prove the description of the simulation data, contrary to 
the early /3-regime, where the microscopic dynamics inter- 
fere and completely mask the critical decay. Similar find- 
ings were reported in |^ for diatomic molecules and in 
||74| for a binary LJ-mixture. 



5 Conclusions 

In this paper we presented results from molecular-dyna- 
mics simulations of a simple model for a supercooled non- 
entangled polymer melt. The temperatures investigated 
are within the range T 0.46 < T < T = 0.52. This range 
is above the critical temperature, Tc = 0.450 ± 0.005, of 
mode-coupling theory (MCT). In order to obtain detailed 
information about the dynamics of the melt we computed 
the incoherent, the collective chain and the collective melt 
intermediate scattering functions as well as their space 
Fourier transforms, the van Hove correlation functions. 
The main focus of this part of the paper was the relax- 
ation of the melt in the /3-regime. The final structural 
a-relaxation is the subject of the subsequent second part 

It was found that the monomer displacements, Vi{t) — 
Yj{t) = 1, . . . ,N), of a chain are not Gaussian dis- 
tributed at all times, as the Rouse model assumes, espe- 
cially when temperature is close to Tc (Figs. |^ and 
However, the Rouse model provides a very good descrip- 
tion of the mean-square displacements {[ri{t) — rj(0)]^). 
This suggests that the Rouse modes of our model remain 
orthogonal at all times, even if the melt is strongly su- 
percooled. The same conclusion was also drawn in 
pm . This property could be important for the further ad- 
vancement of theoretical models for supercooled polymer 
melts (see [|7^, for instance). 

We apphed tests, proposed in the literature fl |65| , |6^ , 
p^,fl], for the space-time factorization theorem of MCT, 
which work directly with the simulation data without in- 
voking any fit procedure (Figs. || and 0). These tests show 
that the theorem is verified in real and reciprocal space at 
all investigated temperatures and that the g-dependence 
of the next-to-leading order correction coefficients quali- 
tatively agrees with MCT predictions. Furthermore, the 
factorization property also seems to hold for the Rouse 
modes. However, higher order corrections behave differ- 
ently than those of the scattering functions (Fig. ||). By 
determining the critical amplitudes for the /3-regime in 
real space conclusions about the typical length scales of 
the dynamics were obtained. The /3-process emerges as a 
localized process which is dominated by the cooperative 
motion of a monomer and its nearest neighbours. The in- 
fluence of the other neighbour shells rapidly decreases to 
zero with increasing distance from the central monomer. 

In a quantitative analysis all coefficients describing the 
correlators in the /3-regime in leading and next-to-lead- 
ing order in the separation parameter were computed and 
their g-dependence was discussed. The analysis suggests 
that there are polymer-specific effects on the length scale 
of the coils, which contribute to the slowing down of struc- 
tural relaxation on the scale of the chains (Figs. || and 
The quantitative comparison with the MCT predictions 
corroborated our qualitative finding that the factoriza- 
tion theorem describes the simulated correlators well in 
the centre of the /3-regime. The parameters of the theo- 
rem, fq and hq, for the coherent scattering functions were 
not fitted, but read off from the simulation data by ap- 
plying general MCT predictions to the results obtained 
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from a p revious study of the incoherent scattering func- 
tion . This shows that MCT can provide guidehnes for 
a quantitative analysis of our simphfied model of a super- 
cooled polymer melt. On the other hand, the next-to-lead- 
ing order corrections have to be fitted. They extend the 
quantitative description of the correlators significantly, es- 
pecially at long times and at large wave- vectors (Fig. p^ , 
but a very precise determination of these correction terms 
over the whole g-range is hard to obtain by fits to the 
simulation data (Fig. RW. 
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